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Abstract 

We consider quarkonium in a iiot QCD plasma which, due to expansion and non-zero viscosity, 
exhibits a local anisotropy in momentum space. At short distances the heavy-quark potential 
is known at tree level from the hard-thermal loop resummed gluon propagator in anisotropic 
perturbative QCD. The potential at long distances is modeled as a QCD string which is screened 
at the same scale as the Coulomb field. At asymptotic separation the potential energy is non- 
zero and inversely proportional to the temperature. We obtain numerical solutions of the three- 
dimensional Schrodinger equation for this potential. We find that quarkonium binding is stronger 
at non-vanishing viscosity and expansion rate, and that the anisotropy leads to polarization of the 
P-wave states. 
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I. INTRODUCTION 



In the era of the Relativistic Heavy Ion Collider at Brookhaven and the Large Hadron 
Collider at CERN the theoretical understanding of in-medium modifications of QCD bound 
states is expected to progress significantly. In this paper we focus on the properties of bound 
states of heavy quarks (quarkonia) in an anisotropic plasma. 

In Quantum Chromodynamics (QCD) with small t' Hooft coupling at short distances 
non- relativistic quarkonium states exist. Their binding energies are much smaller than 
the quark mass rriQ ^ Aqcd {Q = c,b), and their size is much larger than I/itlq. At 
zero temperature, since the velocity of the quarks in the bound state is small, f <^ 1, 
quarkonium can be understood in terms of non-relativistic potential models f]?] using the 
Cornell potential [2| . The potential model can actually be derived directly from QCD as an 
effective field theory (potential non-relativistic QCD - pNRQCD) by integrating out modes 
above the scales mg and then rnqv, respectively 

At high temperatures, the deconfined phase of QCD exhibits screening of static color- 
electric fields [4]. It is expected that this screening leads to the dissociation of quarkonium 
states, which can serve as a signal for the formation of a deconfined quark-gluon plasma 
in heavy ion collisions [5,]. Inspired by the success at zero temperature, potential model 
descriptions have also been applied to understand quarkonium properties at finite temper- 
ature. The pioneering paper of Matsui and Satz pi] was followed by the work of Karsch, 
Mehr and Satz (KMS) |6|, which presented the first quantitative calculation. In recent 
works more involved calculations of quarkonium spectral functions and meson current corre- 
lators obtained from potential models have been performed 0, @, Isl, 10, 11, 12]. The results 
have been compared to first-principle QCD calculations performed numerically on lattices 



13l . |14| . Il5l . |16| . A summary and review of the current understanding of these potential 
models is presented in [l3| , and different aspects of quarkonium in collider experiments can 
be found in [H]. More recently the imaginary part of the potential due to Landau damping 
has been calculated [l^, 20]. Also, the derivation of potential models from QCD via effective 



field theory methods has been extended to finite T [21[. All of these works, however, have 
been performed with the assumption of an isotropic thermal medium. 

Here, we attempt a first assessment of the properties of quarkonium states in a QCD 
plasma which exhibits an anisotropy in momentum space. Such anisotropy may arise due to 
a locally anisotropic hydrodynamic expansion of a plasma with non- vanishing shear viscosity. 
It leads to an angular dependence of the QQ potential ^] . We note that the non-equilibrium 
effect described here arises beyond the linear response approximation in that the operators 
corresponding to various properties of quarkonium states need to be evaluated in an ensemble 
of anisotropic (in momentum space) gauge field configurations. 

After reviewing the anisotropic plasma in Section [Til we formulate the first potential model 
for an anisotropic medium in Section lllli The solution of the 3-dimensional Schrodinger 
equation is described in section [TV] and the numerical results are presented in section |Vl 
We note that by now it is understood that merely solving the Schrodinger equation for the 
individual states is not enough, especially close to the continuum threshold. Here many- 
body interactions should be taken into consideration by solving the Schrodinger equation 
for the non-relativistic Green's function [l2|. However, for the purpose of the present work 
the simple analysis suffices. We summarize our results and draw our conclusions in Section 
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II. THE ANISOTROPIC PLASMA 



The phase-space distribution of gluons is assumed to be given by the following ansatz [2^, 



/(p) = /iso(^Vp2 + e(p-n)2j . (1) 

Thus, /(p) is obtained from an isotropic distribution /iso(|p|) by removing particles with a 
large momentum component along n , the direction of anisotropy. We shall restrict ourselves 
here to a plasma close to equilibrium, which is motivated by the fact that in a heavy-ion 
collision quarkonium states are expected to form when the temperature has dropped to (1- 
2) Tc] by then, the plasma may have equilibrated at least partly. Hence, we assume that the 
function /iso(|p|) is a thermal ideal-gas distribution. 
The parameter ^ determines the degree of anisotropy, 

where p^, = p • n and p^ = p — n(p ■ n) denote the particle momentum along and perpen- 
dicular to the direction n of anisotropy, respectively. If ^ is small then it is also related to 
the shear viscosity of the plasma; for example, for one-dimensional boost-invariant expan- 
sion 
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^^^^ (3) 

Its 

where T is the temperature, r is proper time (and 1/r is the Hubble expansion rate), and 
?7/s is the ratio of shear viscosity to entropy density. In an expanding system, non- vanishing 
viscosity (finite momentum relaxation rate) implies an anisotropy of the particle momenta 
which increases with the expansion rate 1/r. For ~ 0.1 - 0.2 and tT ~ 1 - 3 one finds 
that ^ ~ 1. 

We should stress that in this paper we restrict to solving the time-independent 
Schrodinger equation, i.e. we assume that the plasma is at a constant temperature T and 
anisotropy ^. This approximation is useful if the time scale associated with the bound state, 
~ l/|ii^bind|, is short compared to the time scales over which T and ^ vary. Indeed, for 
sufficiently large quark mass mq this condition should be satisfied. 



III. THE KARSCH-MEHR-SATZ MODEL AT FINITE TEMPERATURE 



A. Isotropic medium = 0) 

Lacking knowledge of the exact heavy-quark potential at finite temperature, different 
phenomenological potentials, as well as lattice-QCD based potentials have been used in 
potential models to study quarkonium. 

The KMS model [6| assumes the following form of the heavy-quark potential at finite 
temperature in an isotropic plasma with ^ = 0: 

F(r, T) = exp (—niD r) H [1 — exp (—rriD r)] . (4) 

r rriD 

Here, a ~ 0.385 is an effective Coulomb coupling at (moderately) short distances, a = 
0.223 GeV^ is the string tension and moiT) is the Debye screening mass. 
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Eq. (jlj) is a model for the action of a Wilson loop of size 1/T and r in the temporal 
and spatial directions, respectively (see [28^] and references therein). This potential has been 
used before to study quarkonium bound states [Z]- However, it was reahzed shortly after 
that eq. (jl]) cannot be taken directly as the heavy-quark potential because it contains an 
entropy contribution; see, for example, the discussion in refs. [l2|, 11, 28, 2^. Rather, eq. 



corresponds to the free energy due to the presence of a QQ in the medium. We emphasize 
that the entropy term in the lattice data is merely a perturbative entropy contribution 
present at large distances^, r — >■ oo, and it is absent at short distances [28|. One can remove 
this entropy term from the lattice data by parameterizing F{r oo,T) = F^oiT) in the 
form Foo(T) = a/T ~ bT and then adding the term bT to F(r, T) at large distance, thereby 



obtaining what has been called the physical potential in [ll|, ll2|. 

Alternatively, one could calculate the full entropy S = —OF/dT and add it to the free 
energy, which leads to the internal energy U = F + TS. The internal energy calculated 
in lattice QCD [30] shows a large increase in Uao near T^, due to the large increase of the 
entropy near T^. Furthermore, at temperatures T ~ a potential model based on the 
internal energy becomes much more binding than the T = Cornell potential (we refer 
to this as "overshooting"). For these reasons, the internal energy U{r,T) obtained on the 
lattice should neither be identified with the heavy quark potential, although it has been 
used in potential models before 0, 10, 2^. Nevertheless, the internal energy provides a 



useful upper limit for the potential at finite T. A version of the internal energy in which the 
overshooting problem was eliminated, was designed in 11, l2| and called the most confining 
potential. 

In this paper we also construct a model for a potential which could be viewed as an upper 
limit for the heavy quark potential, i.e. ~ Uoo- Our present model is very simple and 
contains a minimum number of parameters, as the primary goal is to generalize the finite-T 
potential to anisotropic media. In our model we add the full entropy contribution to the 
KMS ansatz (|4]): 

V{r,T) = Fir,T)-T^^^ (5) 

Q (T 

K, (1 + m^) r) exp {—m£i r) + 2 [1 — exp {—niD r)] — ar exp(— m/) r) . (6) 

r mD 

In the second line we have used that m£, is approximately proportional to T at high tem- 
peratures. Since the effect of the running of the coupling is important only at distances 
< 0.1 fm, not relevant for quarkonium studies, here we do not consider running- coupling 
corrections. Fig. [1] compares the potential at finite temperature to that at mu = which is 
a Cornell potential. 

This potential, just as its original form (jl]), essentially represents an interpolation from 
the well-known Cornell potential at short distance to an exponentially Debye-screened string 
attraction at large r. With (7 ~ 2, mjj — gT and ~ 200 MeV, the length scale where 
medium effects become large is roughly given by rnied(^) — Tc/{2T) fm, in approximate 
agreement with lattice results from ref. jSli- In & corrections to the Cornell potential are 
suppressed at distances r < l/mo, i.e. they appear only at order (ttidt)^. This is due 
to the fact that we subtracted the derivative dF/dlogT even at intermediate distances; it 
appears to give a better representation of the lattice potential at r < rj^^ediT), which in 



We evaluate it in anisotropic HTL resummed perturbation theory in section fill B 21 
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FIG. 1: The model potential from eq. ([6]) at zero and at finite temperature as a function of distance. 
Temperature is normalized to Tc = 192 MeV and the temperature dependence of the Debye mass 
is parameterized as given in eq. (f3T|) below. 



fact coincides with the Cornell ansatz. One can see in Fig. [T] that our potential V{r) is 
very close to the Cornell potential for distances up to r ~ 0.4 fm, in agreement with lattice 
results The finite-temperature potential ([6]) does not overshoot the Cornell potential 
significantly at any r at least up to temperatures on the order of 1.5 Tc. This is actually 



the temperature range where most bound states (except perhaps IS bottomonium) are 



expected to dissociate in an isotropic medium [ll|, ll2| . On the other hand, Fig. [T] shows that 



at rather high temperatures of order 3Tc, the model ([6]) does overshoot the Cornell potential 
at short distances. This indicates that this simple form of the finite-T potential is not 
appropriate when the Debye mass m^i is large. However, this regime is not of interest here 
since even the bb states are no longer bound. Overall, the potential (EI) appears to provide a 
reasonable model for the inter-quark potential in the deconfined phase at (moderately) high 
temperatures. 

At r — > cxD the potential ([6]) approaches 

, , a 0.16 GeV^ 



Again, this is in approximate agreement with the Voo ~ 1/T ansatz used in ref. ll|. Note, 
in particular, that ([7]) is about the same as the internal energy Uoo{T) obtained from the 
lattice data [ll[. We take this as an indication that our potential dHD represents an upper 
limit for the possible finite temperature potentials. 

The main assumption of the KMS model is that the very same screening scale rriD which 
emerges in the Debye-Coulomb potential also appears in the non-perturbative, long-distance 
contribution due to the string. In the following, we take over this assumption to anisotropic 
plasmas. 

It is interesting to note that the KMS ansatz for the free energy from eq. (jl]) can be 
obtained in the usual way from the Fourier transform of the static propagator, provided 
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that a non-perturbative contribution 



(k2 + ml)^ ^ ^ 

is added to the standard HTL resummed propagator [S^]- Here, is a constant of dimen- 
sion two which can be related to the string tension a by matching onto a Cornell potential 
at small m^i r. The presence of such an additional dimensionful scale (besides T) also leads 
to a non- vanishing trace of the energy- momentum tensor (33l |. 



B. Anisotropic medium (^ > 0) 

Our fundamental assumption is that the modified KMS potential ([6]), which provides a 
reasonable upper-limit model for the heavy-quark potential in isotropic media, retains its 
basic form also when the local momentum distribution of the plasma particles is anisotropic 
(^ > 0). However, the isotropic Debye mass m£,{T) is now replaced by an angular dependent 
screening scale /i(6'; ^, T) as discussed in the next section. 



1. Angular dependence of the potential at short distances 



The potential from one-gluon exchange at short distance can be evaluated perturbatively. 
At tree level the potential corresponds to the Fourier transform of the static gluon propagator 
which resums screening effects at high temperature^. For ^ > the potential depends not 
only on the distance between the Q and Q but also on the angle 9 between their separation 
r and the direction n of anisotropy [2^ . 

To linear order in ^, the potential can be expressed as 

/ d^V ,p.r 1 - (p ■ n)Vp2 



nr,0 - ^Ur)-9^MJ (^e-^^-^^-^ (9) 

= VUr){l-mr,e)) , (10) 

where r = rrrio and m£){T) denotes the screening mass in the isotropic medium at a given 

temperature T, as before. Also, V\so{r) is the Debye-screened Coulomb potential in an 
isotropic medium, as given by the first term in eq. (jlj), and the function 

^(f,^^) = /o(f) + /i(f)cos(2^) (11) 

with 



f (^\ _ 6(1 - + f[6 - f(f - 3)] _ f f2 

/o(r) - --6"48+'" ' ^^^^ 



m^ '^-''^^^''^^'^ ^-^^-. (13) 



^ In the first part of tfiis section we are only interested in distances up to ~ It is tfierefore not crucial 
to distinguish carefully whether the potential is identified with the Fourier transform of the propagator 
or with the internal energy. We shall add the entropy contribution later, c.f. eqs. (|18m9p . to suppress 
corrections to the Cornell potential at short distances. 
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Note that eqs. (l9lfT3l) do not apply at large distances f 3> 1, which is a shortcoming of 
the Taylor expansion of the full potential in powers of ^. This is of no importance in the 
following because these expressions are used at short distances only in order to determine 
the angular dependence of the screening scale, which will then replace m£, in eq. 

We can now define the ^-dependent screening mass in the anisotropic medium as the 
inverse of the distance scale rined(^) over which |ry(r)| drops by a factor of e: 

— (X 

log — = 1 , (14) 

/i(^;e,T) = r-Ue;^,T) . (15) 
In f lT^ we have used the fact that rV —a as r — > 0. To leading order in ^ this leads to 

^med = 1 - ^^{fmcd, 0) ■ (16) 

An approximate solution to eq. f|T6l) is given by 

^3 + cos 26* u ^ ^3 + cos26' 

16 m^i 16 

For ^ = 0.5 (1.0) this solution achieves a relative accuracy of < 4% (< 18%) over the entire 
range of 9. The accuracy of this result can be improved systematically by going beyond 
Oif"^) in the expansion of the functions /o and /i introduced in eqs. fll2fl3p . 

The "minimal" extension of the KMS model to non-zero anisotropy consists of replacing 
m£,{T) in ([6]) by fi{9]^,T) from above: 

V{r) = (1 + yur) exp {—fJ.r) H — — [1 — exp {—fir)] — ar exp(— /xr) . (18) 

r /i 

At short distances this reduces to the Cornell potential V{r) = —a/r + ar, as it should be. 

Corrections to eq. (ITS]) due to the finite quark mass can be accounted for by adding a 
temperature- and spin- independent correction proportional to a/{mQr) [34]. This improves 
the accuracy of the wave functions of quarkonium states obtained from the solution of the 
non-relativistic Schrodinger equation. The potential finally takes the form 

V [r) = (1 + yurj exp {—jJ^r) H 1 — exp (— /irj — ar exp(— /ir) ^ — . (19j 

r fx niQ r 

2. Perturbative heavy-quark free energy in an anisotropic medium 

In the limit of infinite mass, the free energy of a heavy quark in a thermal plasma is related 
to the expectation value of Polyakov loops. At high temperature, this can be calculated 
within hard-thermal loop resummed perturbation theory. The leading-order contribution is 
given by 



A-(k)-l 



(20) 



(27r)3 

where A°°(k) is the "hard thermal loop" resummed propagator of static Aq fields. We have 
subtracted the temperature-independent contribution to Fq which is not of interest here. 
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Also, this renders the iritegral UV-finite. For the case of an anisotropic medium, A'"'(k) has 
been calculated in ref. [22|; see, also, ref. [35|. The expression (120|) then turns into 



Fq{^,T) = --a.CpmnhiO , (21) 

with a temperature- independent function h{C,). For small anisotropy the following expansion 
in ^ applies: 

h{^) = 1 - 1^ + l^^e' + ■ • • (22) 
K^J 6^ 240 ^ ^ ^ 

We identify this expression with the free energy of a quark due to its interaction with the 
medium (i.e., half the free energy of a QQ pair at large separation). However, the pertur- 
bative entropy contribution TS = —TdFq/dT should be added again in order to obtain 
the potential. Since (1^ is linear in T (at fixed coupling), it follows that the perturbative 
contribution to Voo vanishes^. On the other hand, lattice data (for an isotropic medium) 
indicate that the free energy of a quark-antiquark pair at infinite separation also contains 
a non-perturbative contribution of the form Vqq{T) = ruied{T) a [l2|], which agrees qualita- 
tively with the prediction of the KMS model, Voo ~ a /mo- Within the framework of the 
KMS model, this implies that at 7^ 0, Voo{0]^,T) depends on angle as screening becomes 
anisotropic. 



IV. SOLVING THE 3D SCHRODINGER EQUATION 

To determine the wave functions of bound quarkonium states, we solve the Schrodinger 
equation 

H = + V{^)+mi+m2 , (23) 

on a three-dimensional lattice in coordinate space with the potential given in eq. f|T9l) . 
Here mi and m2 are the masses of the two heavy quarks and ttir is the reduced mass: 
TTiR = 17111712 / {rrii + 7712) ■ The index v on the eigenfunctions, and energies, E^, represents 
a list of all relevant quantum numbers, e.g. n, I, and m for a radial Coloumb potential. Due 
to the anisotropic screening scale, the wave functions are no longer radially symmetric if 
^ 7^ 0. Since we consider only small anisotropics we nevertheless label the states as IS" 
(ground state) and IP (first excited state), respectively. 

To find solutions to Eq. (123!) we use the finite difference time domain method (FDTD) j36| . 
In this method we start with the time-dependent Schrodinger equation 

z^7/>(x,t) = ^7/>(x,t), (24) 

which can be solved by expanding in terms of the eigenfunctions, 0„: 

V^(x,t) = 5]c„0„(x)e-^^^* . (25) 

V 



^ A term Fq a/T could be generated by a non-perturbative contribution of the form mQ/(k^ + mjj)^ to 
the static gluon propagator, as already mentioned above. rriQ is a constant of dimension two [s^l- 
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If one is only interested in the lowest energy states (ground state and first few excited 
states) an efficient way to proceed is to transform and (l25l) to Euclidean time by a 
Wick rotation, r = it: 

^^(x,r) = -i7V^(x,r), (26) 

and 

7A(x,r) = ^c„0,(x)e-^^^. (27) 

V 



For details of the algorithm we refer to ref . [36[ . 



A. Finding the ground state 

By definition the ground state is the state with the lowest energy eigenvalue, £"0. There- 
fore, at late imaginary time the sum over eigenfunctions fl27p is dominated by the ground 
state eigenfunction 

Jim V'(x,r) ^ Co0o(x)e-^«^ . (28) 

Because of this one can obtain the ground state wavefunction, 0o, and energy, Eq, by solving 
Eq. (l26l) starting from a random three-dimensional wavefunction, ^initial (x, 0), and evolving 
forward in imaginary time. This initial wavefunction should have nonzero overlap with all 
eigenfunctions of the Hamiltonian; however, due to the damping of higher-energy eigenfunc- 
tions at sufficiently late imaginary times we are left with only the ground state, 0o(x). Once 
the ground state wavefunction (or, in fact, any other wavefunction) is found we can compute 
its energy eigenvalue via 

^ „o) = = i^l^ . (29) 

To obtain the binding energy of a state, E^^und, we subtract the quark masses and the 
potential at infinity 

-C^i;,bind = — ITLi — 1712 , , , , , ■ {o\i) 

{(PvWv) 

For the isotropic KMS potential the last term is independent of the quantum numbers v and 
equal to a /mo- In the anisotropic case, however, this is no longer true since the operator 
Voo(^) carries angular dependence, as already discussed above. Its expectation value is of 
course independent of 6 but does depend on the anisotropy parameter ^. 



B. Finding the excited states 

The basic method for finding excited states is to first evolve the initially random wave- 
function to large imaginary times, find the ground state wavefunction, and then project 
this state out from the initial wavefunction and re-evolve the partial-differential equation in 
imaginary time. However, there are (at least) two more efficient ways to accomplish this. 
The first is to record snapshots of the 3d wavefunction at a specified interval Tgnapshot during 
a single evolution in r. After having obtained the ground state wavefunction, one can then 
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go back and extract the excited states by projecting out the ground state wavefunction from 
the recorded snapshots of '?/'(x, r). 

An alternative way to select different excited states is to impose a symmetry condition 
on the initially random wavefunction which cannot be broken by the Hamiltonian evolution. 
For example, one can select the first excited state of the (anisotropic) potential by anti- 
symmetrizing the initial wavefunction around either the x, y, or z axes. In the anisotropic 
case this trick can be used to separate the different polarizations of the first excited state of 
the quarkonium system and to determine their energy eigenvalues with high precision. This 
high precision allows one to more accurately determine the splitting between polarization 
states which are otherwise degenerate in the isotropic Debye-Coulomb potential. 

Whichever method is used, once the wave function of an excited state has been determined 
one can again use the general formulas fl29l) and fl30l) to determine its binding energy. 



V. RESULTS AND DISCUSSION 

In this section we present the solutions of the 3-dimensional Schrodinger equation fl23|) 
in a weakly anisotropic medium. In particular, we determine the temperature dependence 
of the binding energies of different charmonium and bottomonium states obtained with the 
anisotropic potential f[T^ that has been contructed from the most binding isotropic potential. 
The anisotropy- and temperature-dependent screening mass fi{d]C,,T) is given in equation 
f lT7|l . To illustrate the effect of the anisotropy of the medium more clearly we shall also 
compare the results to those obtained for an isotropic medium. In the latter case C = and 
so /i(6';^,T) = mi5(T), where the temperature dependence of the Debye mass is given by 



mDiT)=AgT^{l + nf/Q). (31) 

For Uf = 2 number of massless quark flavors the parameter A = 1.4 has been determined 
in lattice calculations |37|. We choose a fixed gauge coupling of g = 1.72 which yields 
m£){T) /T ^ 2.8. This agrees approximately with lattice estimates of mu/T for temperatures 
on the order of T/TV^ 1.5, and it also gives a reasonable estimate of the free energy at 
infinite separation |37|. The values of the charm and bottom quark masses are chosen such 
that at low temperature the correct masses of Mj/^ = 3.1 GeV and My = 9.4 GeV for the 
J/ip and the T, respectively, are recovered. Accordingly, 

rUc = 1.3 GeV and rub = 4.7 GeV . (32) 

All of the results reported below were obtained from lattices with lattice spacings ap- 
proximately 20 times smaller than the root-mean-square radius (r^)y^(T, ^, mg) of the state 
under consideration, defined by the quantum number v. The lattice size L was chosen to 
be about 6 times larger than the RMS radius'^. Discretization errors and finite-size effects 
are thus expected to be reasonably small and nearly independent of T, ^, mg, and v. We 
stopped the time evolution when the energy Ev{t) of the state had stabilized to within 10~^. 
A more detailed investigation of numerical errors is beyond the scope of the present paper. 



^ Since we restrict the analysis to only weak medium anisotropics, we employ isotropic lattices with uniform 
lattice spacing in all three cartesian directions. 
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FIG. 2: Temperature-dependence of the binding energies l-Ebindl for the ground-states of charmo- 
nium (lower curves) and bottomonium (upper curves) in the vector channel for two values of the 
plasma anisotropy parameter ^. The straight line corresponds to a binding energy equal to the 
temperature. 
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FIG. 3: Expectation value of Voo in the J/^* or T states as a function of temperature for two values 
of the plasma anisotropy parameter ^. 
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Our goal here is to show how quarkonium states may be affected by the anisotropy of the 
medium. 

The temperature dependence of the binding energies of charmonium and bottomonium 
ground states in the vector channel^ are depicted in Fig. [21 The figure shows the results 
obtained for isotropic E_= and anisotropic ^ = 1 media. The former are in agreement with 
those obtained in ref. ll|] with the so-called most- confining isotropic potential. As expected, 
the binding energy decreases as the screening mass m£,{T) increases with temperature T. 
This plot also indicates that |-Ebind| increases with the anisotropy ^. This can be understood 
from the fact that in an anisotropic plasma the screening scale iJ,{6) at a given temperature is 
smaller than the corresponding Debye mass mo] see eq. f|T7j) . As a consequence, the screening 
of the attractive Coulomb and string contributions is less accentuated in the anisotropic 
plasma and so quarkonium states are bound more strongly than in an isotropic medium. 
The magnitude of this effect is substantial even for the moderate anisotropy considered here. 
Near the critical temperature T^, for example, the binding energy of the IS vector cc ground 
state increases by about 50%, and that of the IS bb ground state increases by roughly 30% 
compared to the binding energies calculated in an isotropic medium (the only case addressed 
previously in the literature). 

It is important to highlight another aspect of the reduced screening. In the KMS model 
the asymptotic value of the potential is intrinsically related to the screening mass via the 
relation {Voo)(T) = {v\afj,~^{9; ^,T)\v) . This implies that in the anisotropic medium less 
screening translates into an increase of the potential at infinite separation, Vqo. The above 
is illustrated in Fig. [3] which shows the expectation value of Voo in the T state^. Voo, in turn, 
determines the continuum threshold, which, at a given T is at higher energy than in the 
isotropic case. This implies that at a given temperature the energy gap between the bound 
state and the continuum, which is the binding energy, is increased compared to the isotropic 
case. 

Comparing the behavior of (V"oo) to that of the binding energy of the T from fig. [2] shows 
that the decrease of |£^bind| with T is largely due to the decrease of the continuum threshold 
(Voo). Indeed, we have confirmed that the wave function of the bb ground state is essentially 
unaffected by the medium, i.e. it is almost independent of T (for T ^ 2Tc) and ^ (for ^ ^ 1). 
This has interesting implications for phenomenology: On one hand, the center of the T-peak 
in the dilepton invariant mass distribution may not shift much (since Voo = for decay into 
dileptons) even for temperatures > where the binding energy is significantly lower than 
in vacuum. On the other hand, when [-E'hiriril T we expect substantial broadening of the 



states due to direct thermal activation [111 . |38|. The thermal width can be estimated from 
the binding energy 38|. When the width is larger than the binding energy, a state decays 
faster than it binds ll|. Note, that in the same temperature domain collisions with thermal 
particles of the medium would further broaden the width of a state. Thus, the dissociation 
of the bound states may be expected to occur roughly when |-Ebind| ~ T [17]. With the 
potential investigated in this paper, which likely represents an upper limit for the attractive 
interaction, the condition |£'bind| ~ T is met for the J/ip by 1.2 for ^ = 0, in agreement 



with previous results [llj, and by 1.4Tc for ^ = 1. We stress furthermore that the thermal 



^ Spin effects are neglected in our treatment and the ground state could be identified with either the 
pseudo-scalar or the vector state. For definitness, we shall refer to the vector channel. 

^ We recall that Voo is proportional to the identity at ^ = and hence its expectation value is the same for 
all states. At ^ = f we obtain a very small difference between {J/'i>\Voo\J/'^) and {T\Voc\^)- 
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FIG. 4: Root-mean-square radii of the 15 ground-states of charmonium (upper curves) and bot- 
tomonium (lower curves) as functions of temperature for two values of the plasma anisotropy 
parameter ^. 



density of a given state, 



p ~ exp 



bind 



T 



(33) 



is not enhanced significantly when l-Ebindl < T . In other words, since T decreases with time 
in a heavy-ion collision, quarkonium states with quantum numbers v should appear at a 
temperature T„ ~ |-Ei,,bind| From fig. [2] it is plausible that in a viscous plasma quarkonium 
synthesis occurs at higher temperature than in a perfectly equilibrated medium. For the 
J/^ for example, ATgynth/T'c ~ 20% for ^ = 1 as compared to ^ = 0. 

In fig. m we show the root-mean-square (RMS) radii (r^)^/^(T, ^) of the cc and hb ground 
states as functions of temperature. The former grows rather rapidly about the dissociation 
point where |-Ebind| ~ T. The size of the T, on the other hand, increases only little with 
temperature. We can understand these results, qualitatively, as follows. For charmonium 
the string part of the potential dominates, and the growth of its RMS radius with T indicates 
that screening of the string is strong*. We observe a similar behavior of the first IP excited 
state of bottomonium. On the other hand, IS* bottomonium is too small to be affected 
strongly by screening (for T ^ 2Tc), it is essentially a Coulomb state. The weaker binding 
as compared to low temperature is largely due to a decrease of the continuum threshold Voo-, 
as already mentioned above. 



^ The total number of formed quarkonium states depends on how many heavy quarks are produced in the 
initial hard processes, and on what fraction thereof is bound in D- and _B-mesons, respectively. 

^ Recall that the in the KMS model the string tension enters with a factor of l/mniT) at intermediate 
distances on the order of r ~ l/mD(T). 
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FIG. 5: Absolute energy of the J/^ and T states, minus twice the corresponding quark mass, as 
a function of temperature. 



In this vein, it is also instructive to look at the behavior of the absolute energy of the J/^ 
and T states versus temperature, shown in fig. O We recall that E — 2mQ = (Kx)) + -Ebind- 
The energy of the T increases slightly with temperature as is expected for a small-size state 
bound mainly by the Debye-Coulomb part of the potential (plus a constant): the first term 
on the right hand side of eq. (fT9l) increases with the screening mass fi. On the other hand, 
Ej/^ decreases with T because the second term in (fT9|) decreases as /i increases. 

Fig. [6] shows the temperature dependence of the binding energies of the IP states of 
bottomonium^, identified with the xt- The anisotropy again leads to an increase of |-Ebmd| 
by about 50%, comparable to the behavior of the J/\l/ from above. It also leads to a 
preferred polarization of the Xb, with about 50 MeV splitting between states with angular 
momentum = and = ±1, respectively. At T ~ T^, due to the Boltzmann factor (1331) . 
the population of the state with = is about 30% higher than that of either one of 
the Lz = ±1 states. Here, the polarization is with respect to the axis n of anisotropy, 
which coincides with the direction of expansion of the plasma. In addition, quarkonium 
states produced in high-energy collisions initially through semi-hard gluon fusion may exhibit 
polarization with respect to the particle velocity vector [39| . 



VI. SUMMARY AND CONCLUSIONS 



In a viscous plasma, anisotropic expansion of a fluid element leads to an anisotropy of 
the quasi-particle momentum distributions. The Hard Thermal Loop resummed propagator 



^ The IP states of charmonium do not have binding energies which exceed the temperature significantly 
for T>Tc. 
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FIG. 6: Temperature-dependence of the binding energy for the IP state of bottomonium for two 
values of the plasma anisotropy parameter ^. The straight line corresponds to a binding energy 
equal to the temperature. 



of static chromo-electric fields then carries angular dependence which leads to anisotropic 
Debye screening. 

In this paper we have proposed the first model for the static potential between a very 
heavy quark and anti-quark in a hot anisotropic QCD plasma. Conceptually, we assume 
that the time scale ~ l/|£'bind| associated with the bound state is short compared to the 
time scales over which the temperature and the anisotropy evolve. 

At distances on the order of the Debye length the potential can be calculated from 
perturbative QCD (at high temperature). At larger distances it is however dominated by 
the non-perturbative string attraction. Lattice gauge theory simulations have shown that in 
the deconfined phase the string is screened at a similar scale r-^^^iT) ~ l/m£,{T) and that 
at infinite separation the free energy of a QQ pair approaches a constant V^(T), equal to 
twice the free energy of a single heavy quark in the plasma. 

The essential features appear to be in qualitative agreement with a model proposed by 
Karsch, Mehr and Satz [6] long time ago. However, to obtain the heavy-quark potential we 
subtract the entropy contribution from their ansatz for the free energy of a QQ pair. We 
thereby obtain the internal energy of the QQ pair which should be viewed as an upper limit 
for the physical potential. The latter may be less binding than the KMS internal energy 
used here. We note that in the relevant temperature region, up to about 2Tc, the potential 
at short and intermediate distances follows the zero-temperature Cornell-potential; i.e., the 
overshooting problem of the internal energy is eliminated, in accordance with lattice data 
on the free energy. 

The KMS model correctly reproduces the Cornell potential at short distances and, more- 
over, does not introduce any new parameters besides the string tension. This is important 
for our present goal of extending the isotropic potential to anisotropic plasmas. Knowl- 
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edge of the anisotropic screening scale obtained from the gluon propagator is sufficient to 
generahze the KMS model to anisotropic media. 

We then proceed to solve the Schrodinger equation with this potential to determine the 
wave functions of bound cc and bb states in the plasma. The radial Schrodinger equation is 
no longer sufficient as the potential carries angular dependence. We employ a finite difference 
time domain method (in Euclidean time) on a three-dimensional lattice to obtain the wave 
functions and the binding energies. Some medium effects are neglected in this approach. 
However, solving for the full non-relativistic Green's function (inch threshold effects) in three 
dimensions is beyond the scope of this paper. 

We find that just above the critical temperature Tc ~ 192 MeV for deconfinement (in 
QCD with Nc = 3 colors) in an anisotropic medium both the IS* state of charmonium as 
well as the IS and IP states of bottomonium have binding energies larger than T; the 
temperature may serve as a rough estimate for the width of the states. The binding energies 
decrease with temperature and cease to exceed the estimated width F ~ T at some higher 
temperatures. We note, also, that the Boltzmann enhancement factor exp(— £Jbind/2^) for 
bound states is negligible anyways when |£^bind| ^ T. 

The decrease of |£^bind| with T is due to two effects: First, the continuum threshold 
Voo{T) decreases approximately hke ~ l/T. The energy gap between the bound state and 
the continuum, which is the binding energy, therefore decreases, too. In fact, for the model 
adopted here, we find this to be the dominant effect on the IS ground state of bottomonium 
whose wave function is rather insensitive to the presence of the medium. The state is too 
small to be affected strongly by screening. Hence, the T peak in the dilepton invariant 
mass distribution may not experience a large shift although one should expect substantial 
broadening near the dissociation temperature. 

Larger states such as the 15* ground state of charmonium and the IP excited state of 
bottomonium, however, may also experience some modifications due to screening. The root- 
mean-square radii of these states increase rather rapidly with T around the dissociation 
point l-Ebindl ~ T. 

The two main results of this work are as follows. At fixed T, the screening mass decreases 
with increasing ^. In the KMS model, the asymptotic value of the potential is intrinsically 
related to the screening mass via Voo{0) ~ l/l^'iO; ^, T). Hence, less screening translates into 
an increase of the potential at infinite separation. This implies that the binding energies of 
bound states increase, too. The effect is quite substantial even for moderate anisotropics 
^ ~ 1 considered here: we find that just above Tc the binding energy of the bottomonium 
ground state increases by about 30%, that of 15" charmonium and of IP bottomonium by 
50%. Thus, such quarkonium states may exist up to somewhat higher temperatures than in 
case of an isotropic, perfectly equilibrated medium (for = the J/ip and the T are expected 
to dissociate by 1.2 Tc and 1.8 Tc, respectively, in agreement with previous potential model 
calculations) . 

The other important new effect identified here is that the angular dependence of the 
inter-quark potential in an anisotropic medium induces a polarization of states with non- 
zero angular momentum. According to our estimates, the splitting of the Xb with = 
and Lz = ±1, respectively, is on the order of 50 MeV. At T ~ 200 MeV, the population 
of the state with = is Boltzmann-enhanced by about 30% as compared to the states 
with angular momentum along the direction of anisotropy, respectively. The experimental 
confirmation of such a polarization at RHIC or LHC may provide first evidence for a non-zero 
viscosity of QCD near Tc. 
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The next step of our investigation is the determination of the imaginary part of the 
potential, which will provide insight into how the anisotropy of the medium affects the 
widths of the states. We shall present results in a future publication. 
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